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The quark-model baryon-baryon interaction fss2, proposed by the Kyoto-Niigata group, 
is a unified model for the complete baryon octet (Bg = N, A, £ and S), which is formu- 
lated in a framework of the (3q)-(3q) resonating- group method (RGM) using the spin-flavor 
SUe quark-model wave functions and effective meson-exchange potentials at the quark level. 
Model parameters are determined to reproduce properties of the nucleon-nucleon system 
f*^ ' and the low-energy cross section data for the hyperon-nucleon scattering. Due to the sev- 

| eral improvements including the introduction of vector-meson exchange potentials, fss2 has 

, achieved very accurate description of the AW and YN interactions, comparable to various 

■ one-boson exchange potentials. We review the essential features of fss2 and our previous 

model FSS, and their predictions to few-body systems in confrontation with the available 
| experimental data. Some characteristic features of the BgB$ interactions with the higher 

, strangeness, S = —2, — 3, — 4, predicted by fss2 are discussed. These quark-model inter- 

actions are now applied to realistic calculations of few-body systems in a new three-cluster 
Faddeev formalism which uses two-cluster RGM kernels. As for the few-body systems, we 
discuss the three-nucleon bound states, the ANN-SNN system for the hypertriton, the aaA 
system for ^Be, and the AAa system for yf^He. 
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§1. Introduction 



^ ■ An important purpose of studying baryon-baryon interactions in the quark model 

(QM) is to obtain the most accurate understanding of the fundamental strong in- 
teractions in a natural picture, in which the short-range part of the interaction is 
relevantly described by the quark-gluon degrees of freedom and the medium- and 
long-range parts by dominated meson-exchange processes. In the spin-flavor SUq 
QM, the baryon-baryon interactions for all the octet baryons (Bq = N, A, U and 
S) are consistently treated with the well-known nucleon-nucleon (NN) interaction. 
We have recently proposed a comprehensive QM description of general baryon-octet 
baryon-octet (B$Bs) interactions, which is formulated in the (3q)-(3q) resonating- 
group method (RGM) using the spin-flavor SUq QM wave functions, a colored version 
of the one-gluon exchange Fermi-Breit interaction, and effective meson-exchange po- 
tentials (EMEP's) acting between quarks. 1 )" 5 ) The early version, the model FSS, 1 )" 3 ) 
includes only the scalar (S) and pseudoscalar (PS) meson-exchange potentials as the 
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EMEP's, while the renovated one fss2 4 )' 5 ) introduces also the vector (V) meson- 
exchange potentials and the momentum-dependent Bryan-Scott terms for the S and 
V mesons. Owing to these improvements, the model fss2 in the NN sector has at- 
tained the accuracy comparable to that of one-boson exchange potentials (OBEP's). 

These QM potentials can now be used for various types of many-body calcula- 
tions, which include the G-matrix calculations 6 ) of baryonic matter and the Faddeev 
calculations of few-baryon systems. 7 -* For this purpose, we have recently developed 
a new three-cluster formalism which uses two-cluster RGM kernels explicitly. 8 ) The 
proposed equation eliminates three-cluster redundant components by the orthogonal- 
ity of the total wave function to the pairwise two-cluster Pauli-forbidden states. The 
explicit energy dependence inherent in the exchange RGM kernel is self-consistently 
determined. This equation is entirely equivalent to the Faddeev equation which 
uses a modified singularity-free T-matrix (which we call the RGM T-matrix) gener- 
ated from the two-cluster RGM kernel. We first applied this formalism to a three- 
dineutron system and the 3a system, and obtained a complete agreement between 
the Faddeev calculations and variational calculations which use the translationally 
invariant harmonic-oscillator (h.o.) basis. 8 )' 9 ) Here we apply the formalism to the 
Faddeev calculations of the three-nucleon (3N) bound state, 7 ) and the ANN-UNN 
system for the hypertriton, as well as the aaA and AAa systems. 

§2. Formulation 

The present model is a low-energy effective model which introduces some of the 
essential features of QCD characteristics. The color degree of freedom of quarks is 
explicitly incorporated into the non-relativistic spin-flavor SUq quark model, and 
the full antisymmetrization of quarks is carried out in the RGM formalism. The 
gluon exchange effect is represented in the form of the quark-quark interaction. The 
confinement potential is a phenomenological r 2 -type potential, which has a favor- 
able feature that it does not contribute to the baryon-baryon interactions. We use 
a color analogue of the Fermi-Breit (FB) interaction, motivated from the dominant 
one-gluon exchange process in conjunction with the asymptotic freedom of QCD. 
We postulate that the short-range part of the baryon-baryon interaction is well de- 
scribed by the quark degree of freedom. This includes the short-range repulsion and 
the spin-orbit force, both of which are successfully described by the FB interaction. 
On the other hand, the medium range attraction and the long-range tensor force, 
especially afforded by the pions, are extremely non-perturbative from the viewpoint 
of QCD. These are therefore most relevantly described by the effective meson ex- 
change potentials (EMEP). In the previous model called FSS, 1 )" 3 ) only the scalar 
and pseudoscalar mesons are introduced. The model fss2 4 )' 5 ) is the most advanced 
model which also includes the vector-meson exchange EMEP. The full QM Hamil- 
tonian H consists of the non-relativistic kinetic-energy term, the phenomenological 
confinement potential Ujj , the colored version of the full FB interaction Uf^ with ex- 
plicit quark-mass dependence, and the EMEP U[j generated from the scalar (4?=S), 
pseudoscalar (PS) and vector (V) meson exchange potentials acting between quarks: 
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H = t(^ 2 + £-.- T o)+t(uS' + ^ + u! 3 + U^ + U^. (2-1) 

i=i v » / j <Jr - 

The RGM equation for the relative-motion wave function x( r ) reads 

(4>(3q)<P(3q)\E - H\A{<t>(3q)</>(3q) X (r)}} = 0. (2-2) 

We solve this RGM equation in the momentum representation. 10 ) If we rewrite the 
RGM equation in the form of the Schrodinger-type equation as [e — Hq — Vrgm(e)] 
x x( r ) = 0) the potential term, Vrgm(e) = Vd+G+sK, becomes nonlocal and energy 
dependent. Here Vb represents the direct potential of EMEP's, G includes all the 
exchange kernels for the interaction and kinetic-energy terms, and K is the exchange 
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Fig. 2. Same as Fig. 1 but for the higher partial waves with J < 4. 

normalization kernel. We calculate the plane- wave matrix elements of Vr.gm(£)> and 
set up with the Lippmann-Schwinger equation of the RGM T-matrix. This approach 
is convenient to proceed to the G-matrix calculations. 6 )' n ) Faddeev calculations 
using these exchange kernels are also possible with some special considerations of 
the Pauli forbidden states. 8 )' 9 ) 



§3. B 8 B 8 interactions by fss2 and FSS 

3.1. NN properties and G-matrix calculations of nuclear matter 

Figures 1 and 2 display some important low-partial wave NN phase-shift pa- 
rameters predicted by the model fss2, in comparison with the phase shift analysis 
SP99. 12 ) The previous results by FSS are also shown with the dotted curves. Due 
to the inclusion of V mesons, the NN phase shifts of the fss2 at the non-relativistic 
energies up to Ti a b = 350 MeV are greatly improved, and now have attained the 
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Fig. 3. Nuclear matter saturation curves ob- 
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sults of the Paris potential 17 - 1 and the Bonn 
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the non-relativistic calculation in Ref. 19). 
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Fig. 4. Same as Fig. 3 but for the neutron mat- 
ter saturation curve predicted by fss2. The 
normal density po corresponds to the Fermi 
momentum k F — 1.35 fm -1 . The results of 
the Nijmegen soft-core potential (NSC97) 20) 
and the Bonn-B potential 18 - 1 are also shown 
for comparison. 



accuracy almost comparable to that of OBEP's. The good reproduction of the NN 
phase-shift parameters in fss2 continues up to T| a b ~ 600 MeV, 4 ) where the inelas- 
ticity of the 5-matrix becomes appreciable. 

For the correct evaluation of the triton binding energy, it is well known that the 
proper reproduction of the deuteron D-state probability and the 1 So effective range 

Table I. The deuteron properties by fss2 and FSS in the isospin basis. The results by the Bonn-B 
potential 18 - 1 are also shown for comparison. A small difference in FSS from Table IV of Ref. 2) 
is due to the numerical inaccuracy in the previous calculation. The effect of the meson exchange 



current is not included in 


the calculations of Q d and jid- 








FSS 


fss2 


Bonn B 


Expt. 


Ref. 


e d (MeV) 


2.2561 


2.2247 


2.2246 


2.224644 ± 0.000046 


13) 


Pd (%) 


5.86 


5.49 


4.99 






r) = A D /A s 


0.0267 


0.0253 


0.0264 


0.0256 ± 0.0004 


14) 


rms (fm) 


1.963 


1.960 


1.968 


1.9635 ±0.0046 


13) 


Qd (fm 2 ) 


0.283 


0.270 


0.278 


0.2860 ±0.0015 


15) 


Hd (mjv) 


0.8464 


0.8485 


0.8514 


0.857406 ± 0.000001 


16) 
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parameters is essential. Table I shows the deuteron properties predicted by FSS and 
fss2, in comparison with the experiment. The predictions by the Bonn-B potential, 
which has a smaller deuteron .D-state probability, is also shown for comparison. 

Figure 3 shows saturation curves of symmetric nuclear matter, obtained by fss2 
and FSS. 6 ) They depend on the prescription how one deals with the energy spectrum 
of the intermediate single-particle (s.p.) states. In both cases of the QTQ and 
the continuous choices, they are very similar to the predictions by meson exchange 
potentials. In the continuous choice, the saturation curves of fss2 and FSS are very 
similar to that of the Bonn-B potential. 18 ) Unfortunately, our results share the 
common unsatisfactory feature of any non-relativistic models, that the saturation 
point does not deviate from the Caester line. 

The saturation curve for neutron matter, predicted by fss2, is also shown in 
Fig. 4. The results by the Nijmegen soft-core potential (NSC97) 20 ) and the Bonn-B 
potential 18 ^ are also shown for comparison. Since the strongly attractive 3 Si + 3 D\ 
channel does not exist in this case, the total energy per neutron in neutron matter 
becomes repulsive in any models. We find that our quark model potential gives 
very similar results to the standard meson-exchange potentials, as long as the NN 
interaction is concerned. 

The s.p. potentials of the nucleon and hyperons in symmetric nuclear matter are 
shown in Figs. 5 (&f = 1.35 fm -1 ) and 6 (fcjr = 1.07 fm -1 ) for the model fss2. For 
the standard Fermi momentum /cp = 1.35 fm -1 , which corresponds to the normal 
density po, these are fairly deep potentials. For the comparison with the depth of 
the s.p. potentials in the finite nuclei, one has to take a smaller value for fcp because 
of the surface effects. If we assume kF = 1-07 fm" 1 , the s.p. potentials become much 
shallower, and the depth becomes almost comparable with the empirical values, —50 
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Fig. 7. Calculated YW total cross sections compared with the available experimental data. 

MeV for N, -30 MeV for A, and -10 14 MeV for ~. The sign of the E s.p. 

potential is still controversial. 

3.2. YN and YY interactions by fss2 and FSS 

The total cross sections of the hyperon-nucleon (YN) scattering predicted by 
fss2 are compared with the available experimental data in Fig. 7. The "total" cross 
sections for the scattering of charged particles (i.e., E + p and E~p systems) are cal- 
culated by integrating the differential cross sections over cos0 m i n = 0.5 ~ cos6* max = 
—0.5. The solid curves indicate the result in the particle basis, while the dashed 
curves in the isospin basis. In the latter case, the effects of the charge symmetry 
breaking, such as the Coulomb effect and the small difference of the threshold en- 
ergies for E~p and E°n channels, are neglected. New experimental data for E~p 
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elastic total cross sections at the intermediate energies, ps = 400-700 MeV/c, mea- 
sured at KEK, 21 ) are consistent with the fss2 predictions. 

The U + p differential cross sections at the intermediate energy pz = 450 MeV/c 
are compared with the KEK experiment 22 ) in the left panel of Fig. 8. We need more 
experimental data to increase the statistics, in order to see which model is the most 
appropriate. In the right panel, the polarization observables for the U + p elastic 
scattering at pz = 800 MeV/c are shown for the models fss2 and FSS. The recent 
experimental data from KEK-PS 45 7 23 ) imply the asymmetry parameter a exp = 
0.44 ±0.2 at pe = 800 ±200 MeV/c, which is not inconsistent with our quark-model 
predictions although the specific scattering angle is not possible to be identified. 

We show some comparison of the E~ N total cross sections which are recently ob- 
tained from the BNL-E906 experiment. 24 ) The in-medium E~ N total cross sections 
around the momentum region pi^ ~ 550 MeV/c are estimated as cr=Ar(in medium) = 
30 ± 6.7^3'g mb. More detailed analysis using the Eikonal approximation 25 ) gives 
asN(in medium) = 20.9 ± 4.5^2 4 m b- They have also estimated the cross section 
ratio cr~- p /&s-n = l-l^o 7-0 4- ^ we compare these experimental data with the FSS 
and fss2 predictions depicted in Fig. 9, we find that the FSS predictions in the left 
panels seem to be more favorable. However, we definitely need more experimental 
data with higher statistics. 

3.3. Characteristics of the BqBs Interactions by fss2 

Since our quark-model parameters are fixed by using the experimental data in 
the iViV and YN sectors, the B$B$ interactions beyond the strangeness S = —1 
are all model predictions. 5 ) For the systematic understanding of these interactions, 
it is convenient to discuss them by using the SU3 representation basis for the two- 
baryon systems. This is because the quark-model Hamiltonian is approximately 
SU3 scalar, and the interactions with the same SU3 label (A//) should have very 
similar characteristics, as long as the flavor symmetry breaking is negligible. Table II 
illustrates how the two baryon systems in the isospin basis are classified as the 
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Fig. 9. S~p and S°p total cross sections predicted by FSS (left) and fss2 (right). For E~p cross 
sections both of the isospin 1 = and 1 channels contribute, while for E~n (or S°p) only the 
channel with 7 = 1 contributes. The Coulomb force is neglected. 

superposition of the SU3 basis. It has entirely different structure between the flavor 
symmetric and antisymmetric cases. In the 1 Sq state, for example, there appear 
many states having the dominant (22) components. The S-state BgB$ interactions 
for these states should be very similar to that for the NN 1 So state. The (H) s 
component is completely Pauli forbidden and is characterized by the strong repulsion 
originating from the quark Pauli principle. The (00) component in the H-particle 
channel is attractive from the color-magnetic interaction. On the other hand, in the 
NN sector the 3 Si state with I = is composed of the pure (03) state, and the 
deuteron is bound in this channel owing to the strong one-pion tensor force. This 
SU3 state in the flavor antisymmetric case is converted to the (30) state in the larger 
side of the strangeness. Since the (30) state is almost Pauli forbidden, the interaction 
is strongly repulsive. Therefore, the EE interaction is not so attractive as NN, since 
they are combinations of (22) and (30). The other SU3 state (ll) a turns out to have 
very weak interaction. After all, the strangeness S = — 2 sector is most difficult, 
since it is a turning point of the strangeness. It is also interesting to see that EE 
channel with 1 = 3/2 should be fairly attractive, since the same (22) and (03) SU3 
states as in the NN system appear in this isospin channel. 

Figure 10 shows fss2 predictions of the 1 Sq phase shifts for various B$B$ inter- 
actions having the pure (22) configuration. Although the UU interaction with the 
isospin / = 2 is very similar to the NN interaction, the other interactions generally 
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get weaker as the strangeness involved becomes larger. This is a combined effect of 
the flavor symmetry breaking in the quark and meson-exchange contributions. In 
particular, the EE interaction has the lowest rise of the phase shift, which is less 
than 30°. Accordingly, the the EE total cross sections become much smaller than 
the other systems. We find that the NN interaction is the strongest and has the 
largest cross sections among any combinations of the octet baryons. 

Figure 11 shows the 1 So phase shift curves, predicted by fss2, in the AA-EN- SU 
coupled-channel system with the isospin 1 = 0. The maximum peak of the A A phase 
shift is less than 20°, which is much smaller than the previous predictions by various 
models. This result is in good agreement with the recent experimental data for the 
double A hypernucleus yf^He. The finding of this event, called the Nagara event, 26 ) 
is one of the most important contributions in recent years, since the assignment 
of the decaying scheme is very definite. The AB AA value, defined by AB AA = 
B AA ( A 6 A Jie) - 2B A (%H.e), is deduced from this event as AB AA = 1.01 ± 0.20 MeV, 
which implies a weak attraction for the AA interaction. Our G-matrix calculation of 
this system yields an almost right answer, AB AA = 1.12 ~ 1.24 MeV, by taking into 
account the important contribution from the a-particle rearrangement energy. 27 ) In 

Table II. The relationship between the isospin basis and the flavor-S^ basis for the BgBg, systems. 



The flavor- SU3 symmetry is given by the Elliott notation (A/i). The heading V denotes the flavor 
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the next section, we will solve the Faddeev equation of the AAa system, by using 
the AA T-matrix predicted by our QM interaction. 

Some of the following characteristics of fss2 for the BgBg interactions are very 
much different from the Nijmegen predictions given by Stoks and Rijken. 28 ) 



o There is no bound state in the BgBg system, except for the deuteron. 
o The EE total cross sections are not so large as the NN total cross sections, 
o The EN interaction has a strong isospin dependence like the EN interaction, 
o The E~E~ (namely, EE(I = 3/2)) interaction is fairly attractive. 



§4. Faddeev calculation 



4.1. Three-nucleon bound state 

Since our QM B$Bg interaction describes the short-range repulsion very differ- 
ently from the meson-exchange potentials, it is interesting to examine the three- 
nucleon system predicted by fss2 and FSS. Here we solve the Faddeev equation for 
3 H, by directly using the QM RGM kernel in the isospin basis. 7 ) Our Faddeev cal- 
culation is the full 50-channel calculation up to the maximum angular momentum 
J = 6, and the values are almost completely converged as seen in Table III. 29 ) The 
fss2 prediction, —8.52 MeV, seems to be overbound in comparison with the exper- 
imental value S exp ( 3 H) = —8.48 MeV. In fact, this is not the case, since all these 
calculations neglect the charge dependence of the NN interaction, in which the 1 So 
interaction of the nn system is less attractive than that of the np system. The ef- 
fect of the charge dependence is estimated to be —0.19 MeV for the triton binding 
energy. 18 ) If we take this into account, our result is still 150 keV less bound. If we at- 
tribute this difference to the effect of the three-nucleon force, it is by far smaller than 
the generally accepted values, 0.5 ~ 1 MeV, 30 ) predicted by many Faddeev calcu- 
lations employing modern realistic meson-theoretical NN interactions. The charge 
root-mean-square radii of 3 H and 3 He are correctly reproduced. An important point 
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Table III. The three-nucleon bound state properties predicted by the Faddeev calculation with fss2 
and FSS. The np interaction is used in the isospin basis. The heading "No. of channels" implies 
the number of two-nucleon channels included; n max is the dimension of the diagonalization for the 
Faddeev equation; ejvjv is the NN expectation value determined self-consistently; E( H) is the 
ground-state energy; and \J (r 2 )3 H (\J {r 2 )3 Hc ) is the charge rms radius for 3 H ( 3 He), including 
the proton and neutron size corrections. The Coulomb force and the relativistic corrections are 
neglected. 



model 


No. of 


^max 


ejvjv 


E( 3 B) 


VV>3h 






channels 




(MeV) 


(MeV) 


(fm) 


(fm) 




2 (S) 


2,100 


2.361 


-7.807 


1.80 


1.96 




5 (SD) 


5,250 


4.341 


-8.189 


1.75 


1.92 


fss2 


10 (J < 1) 


10,500 


4.249 


-8.017 


1.76 


1.94 




18 (J < 2) 


18,900 


4.460 


-8.439 


1.72 


1.90 




34 (J < 4) 


35,700 


4.488 


-8.514 


1.72 


1.90 




50 (J < 6) 


112,500 


4.492 


-8.519 


1.72 


1.90 




2 (S) 


2,100 


2.038 


-7.675 


1.83 


1.99 




5 (SD) 


5,250 


3.999 


-8.034 


1.78 


1.95 


FSS 


10 (J < 1) 


10,500 


3.934 


-7.909 


1.78 


1.97 




18 (J < 2) 


18,900 


4.160 


-8.342 


1.74 


1.93 




34 (J < 4) 


35,700 


4.175 


-8.390 


1.74 


1.92 




50 (J < 6) 


112,500 


4.177 


-8.394 


1.74 


1.92 



here is that we can reproduce enough binding energy of the triton, without reducing 
the deuteron D-state probability. 

The self-consistent energy of the two-cluster RGM kernel, £nn, has a clear 
physical meaning related to the decomposition of the total triton energy E into the 
kinetic-energy and potential-energy contributions: enn = E/3 + (Ho)/ 6. Table IV 
shows this decomposition, together with the results of CD-Bonn and AV18 poten- 
tials. 30 ) We find that our quark model results by FSS and fss2 are just middle 
between these potentials, which have very different strengths of the tensor force. 

4.2. The hypertriton 

Next, we apply our QM NN and YN interactions to the hypertriton QH) with 
the small separation energy of the yl-particle, Ba bxp = 130 ± 50 keV. 31 ) Since the 
yl-particle is far apart from the two-nucleon subsystem, the on-shell properties of 

Table IV. Decomposition of the total triton energy E into the kinetic-energy and potential-energy 
contributions: E — (Ho) + (V). The unit is in MeV. In the present framework, this is given by 
the expectation value ejvjv of the two-cluster Hamiltonian with respect to the Faddeev solution, 
which is determined self-consistently. The results of CD-Bonn and AV18 are taken from Ref. 30). 



model 


£nn 


E 


(Ho) 


(V) 


fss2 


4.492 


-8.519 


43.99 


-52.51 


FSS 


4.177 


-8.394 


41.85 


-50.25 


CD-Bonn 


3.566 


-8.012 


37.42 


-45.43 


AV18 


5.247 


-7.623 


46.73 


-54.35 
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Table V. Results of the hypertriton Faddeev calculations by fss2 and FSS. The heading E is the 
hypertriton energy measured from the ANN threshold; Ba is the A separation energy; ejvjv 
(ban) is the AW (AN) expectation value determined self-consistently; and Pe is the SNN 
probability in percent. The norm of admixed redundant components is less than 10~ 9 . 



model 


No. of 


E 


B A 


Sn n 


EAN 


Pe 




channels 


(MeV) 


(keV) 


(MeV) 


(MeV) 


(%) 




6 (S) 


-2.362 


137 


-1.815 


3.548 


0.450 




15 (SD) 


-2.423 


198 


-1.762 


5.729 


0.652 


fss2 


30 (J < 1) 


-2.403 


178 


-1.786 


5.664 


0.615 




54 (J < 2) 


-2.498 


273 


-1.673 


5.974 


0.777 




102 (J < 4) 


-2.513 


288 


-1.658 


6.022 


0.804 




150 (J < 6) 


-2.514 


289 


-1.657 


6.024 


0.805 




6(5) 


-2.910 


653 


-1.309 


3.984 


1.022 




15 (SD) 


-2.967 


710 


-1.433 


6.171 


1.200 


FSS 


30 (J < 1) 


-2.947 


691 


-1.427 


6.143 


1.191 




54 (J < 2) 


-3.121 


865 


-1.323 


6.467 


1.348 




102 (J < 4) 


-3.134 


877 


-1.317 


6.488 


1.360 




150 (J < 6) 


-3.134 


878 


-1.317 


6.488 


1.361 



the AN and UN interactions are expected to be well reflected in this system. In 
particular, this system is very useful to learn the relative strength of 1 5o and 3 Si 
attractions of the AN interaction, since the component plays a more important 
role than 3 Si in this system and the available low-energy Ap total cross section data 
cannot discriminate many possible combinations of the 1 Sq and 3 Si interactions. 
In fact, Ref. 32)-34) showed that most meson-theoretical interactions fail to bind 
the hypertriton, except for the Nijmegen soft-core potentials NSC89, 35 ) NSC97f and 
NSC97e. 20 ) It is also pointed out in Refs. 32) and 34) that a small admixture of the 
SNN components less than 1% is very important for this binding. We therefore 
carry out the ANN-ENN coupled-channel Faddeev calculation, by properly taking 
into account the existence of the SU3 Pauli forbidden state (ll) s at the quark level. 

Table V shows the results of the Faddeev calculations using fss2 and our previous 
model FSS. In the 15-channel calculation including the S and D waves of the NN 
and YN interactions, we have already obtained Ba = — £d ~ E( A H.) ~ 200 keV for 
fss2. The convergence with the extension to the higher partial waves is very rapid, 
and the total angular-momentum of the baryon pairs with J < 4 is good enough for 
1 keV accuracy. As for the converged Ba values with 150-channel ANN and ENN 
configurations, we obtain Ba = 289 keV and the ENN component Pjj = 0.80% for 
the fss2 prediction, and B A = 878 keV and P E = 1.36% for FSS. 

Table VI shows the correlation between the A separation energy Ba and the 
1 Sq and 3 S*i effective range parameters of FSS, fss2 and NSC89 AN interactions. 
Although all of these AN interactions reproduce the low-energy AN total cross 
section data below p A ~ 300 MeV/c within the experimental error bars, our quark- 
model interactions seem to be slightly more attractive than the Nijmegen soft-core 
potential NSC89. 35 ) The model FSS gives a large overbinding since the 1 S'o AN 
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Table VI. 1 5'o and 3 Si effective range parameters of various AN interactions and the A separation 
energy Ba of the hypertriton. The Ba value for NSC89 is taken from Ref. 34). 



model 


a s (fm) 


r s (fm) 


a t (fm) 


n (fm) 


B A (keV) 


FSS 2) 


-5.41 


2.26 


-1.03 


4.20 


878 


fss2 4 >' 5 > 


-2.59 


2.83 


-1.60 


3.01 


289 


NSC89 35) 


-2.78 


2.88 


-1.41 


3.11 


143 



interaction is strongly attractive. The phase-shift difference of the 1 So and 3 Si 
states at p A ~ 200 MeV/c is S^So) - <5( 3 Si) ~ 29° for FSS, while ~ 7° for fss2. 
Since the present fss2 result is still slightly overbound, this difference should be 
somewhat smaller in order to reproduce the correct experimental value E>A exp = 
130 ± 50 keV. From the two results given by fss2 and FSS, we can extrapolate the 
desired difference to be 0° ~ 2°, which is consistent with the result in Ref. 36) using 
simulated interactions of the Nijmegen models. 

Table V also shows that the expectation value of the NN Hamiltonian, snn, 
determined self-consistently is rather close to the deuteron energy —£d, especially 
in fss2. This feature is even marked if we decompose these energies to the kinetic- 
energy and potential-energy contributions. Table VII shows this decomposition with 
respect to fss2, FSS and NSC89. (For this comparison, we use the definition of the 
kinetic-energy part of the deuteron, given by h d = {xd\fiNN\Xd)/(Xd\Xd), where Xd 
is the RGM relative wave function between the neutron and the proton.) In fss2, the 
kinetic-energy of the NN subsystem is 1.88 MeV larger than that of the deuteron, 
which implies that the NN subsystem shrinks by the effect of the outer /1-particle, 
in comparison with the deuteron in the free space. In NSC89, this difference is even 
smaller; i.e., 1.18 MeV. These results are consistent with the fact that the hypertriton 
in NSC89 is more loosely bound (B A = 143 keV) 34 ) than in fss2 (289 keV), and 
the yl-particle is very far apart from the NN cluster. The UNN probability in 
NSC89 is P s = 0.5 %. 32 ) Table VII also lists the kinetic -energy and potential- 
energy decompositions for the AN -EN averaged YN expectation value eyn and 
the total energy E. The kinetic energies of eyn are much smaller than those of 



Table VII. Decomposition of the NN and YN expectation values (enn and £yn), the deuteron 
energy (— £<j) and the total three-body energy E to the kinetic-energy and potential-energy 



contributions. 


The unit is ii 


i MeV. 


The results for NSC89 


are taken 


from Ref 


.32). 


model 


tlNN 


+ Vnn 


= £jvjv 


h d + Yd = 


= — Ed (deuteron) 


FSS 


19.986 - 


- 21.303 


= -1.317 


16.982 - 


19.238 = 


-2.256 


fss2 


19.376 - 


■ 21.032 


= -1.657 


17.495 - 


19.720 = 


-2.225 


NSC89 


20.48 - 


- 22.25 


= -1.77 


19.304 - 


21.528 = 


-2.224 


model 


tlYN 


+ Vyn 


= £yn 


(H) + (V) = 


E 


FSS 


10.036 


- 4.602 


= 5.435 


27.372 - 


30.506 = 


-3.134 


fss2 


8.071 


-2.671 


= 5.401 


23.860 - 


26.374 = 


-2.514 


NSC89 


7.44 


- 3.54 = 


= 3.90 


23.45 - 


• 25.79 = 


-2.34 
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&NN i which indicates that the relative wave functions between the hyperon and the 
nucleon are widely spread in the configuration space. The comparison of the total- 
energy decomposition shows that the wave functions of fss2 and NSC89 may be very 
similar. A clear difference between fss2 and NSC89 appears in the roles of higher 
partial waves. The energy increase due to the higher partial waves than the S and 
D waves is 91 keV in fss2 and 168 keV in FSS, respectively. On the other hand, the 
results in Refs. 32) and 34) imply that this is only 20 - 30 keV in the case of NSC89. 
This difference can originate from both of the NN and YN interactions. Since the 
characteristics of the meson-theoretical YN interactions in higher partial waves are 
a priori unknown, more detailed analysis of the fss2 results might shed light on the 
adequacy of the quark-model baryon-baryon interactions. 

4.3. The aaA system for \Be 

As another typical example of three-cluster systems composed of two identical 
clusters, we apply the present formalism to the aaA Faddeev calculation for 9 Be, 
using the aa RGM kernel and the Aa folding potential generated from a simple 
AN effective interaction. 37 ) For the aa RGM kernel, we use the the three-range 
Minnesota force 38 ) with the Majorana exchange mixture u = 0.94687, and the h.o. 
width parameter, v = 0.257 fm~ 2 , assumed for the (0s) 4 a-clusters. The aa phase 
shifts are nicely reproduced in the aa RGM calculation, using this effective NN 
interaction. For the 3a system, we find that the ground-state energy obtained by 
solving the present 3a Faddeev equation is only 1.8 MeV higher than that of the 
fully microscopic 3a RGM calculation. The effective AN interaction, denoted by 
SB (Sparenberg-Baye potential) in Table VIII, is constructed from the 1 S'o and 3 Si 
phase shifts predicted by the YN sector of the model fss2, by using an inversion 
method based on supersymmetric quantum mechanics. 39 -* These are simple two- 
range Gaussian potentials which reproduce the low-energy behavior of the AN phase 
shifts in Fig. 12, obtained in the full coupled-channel calculations: 

Vi So (r) = -128.0 exp(-0.8908 r 2 ) + 1015 exp(-5.383 r 2 ) , 
V 3Sl (r) = -56.31/ exp(-0.7517 r 2 ) + 1072 exp(-13.74 r 2 ) , 
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■a 
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-60 
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800 1000 



Fig. 12. AN- EN 1 S (a) and 3 Si (b) phase shifts for the isospin I = 1/2 channel, calculated with 
fss2 (solid and dashed curves) and with the SB potential (circles). 
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Table VIII. The ground-state energy E gT (0 + ) and 2 + excitation energy E^(2 + ) in MeV, calculated 
by solving the Faddeev equation for the aaA system in the LS coupling scheme. The aa RGM 
kernel is generated from the three-range Minnesota force. The AN force, SB, stands for the 
Sparenberg-Baye potential and NS - JB are the G-matrix based effective forces used by Hiyama 
et al. 42 ^ In Ref. 42), the orthogonality condition model (OCM) is used for the aa interaction. 



The experimental values 


are cited from Ref. 43). 




Van 


£ gr (0+) (MeV) 


£ x (2+) (MeV) 




present 


Ref. 42) 




SB 


-6.837 




2.915 


NS 


-6.742 


-6.81 


2.916 


ND 


-7.483 


-7.57 


2.935 


NF 


-6.906 


-7.00 


2.930 


JA 


-6.677 


-6.76 


2.919 


JB 


-6.474 


-6.55 


2.911 


Exp't 




-6.62 ±0.04 


3.029(3) 








3.060(3) 



where V(r) in MeV and r in fm is the relative distance between A and N. In the 
3 Si state, the phase-shift behavior only in the low-momentum region with p\ a ^ < 
600 MeV/c is fitted, since the cusp structure is never reproduced in the single- 
channel calculation. Since any central and single-channel effective AN force leads 
to the well-known overbinding problem of ^He by about 2 MeV (in the present 
case, it is 1.63 MeV), 40 ) the attractive part of the 3 S± AN potential is modified 
to reproduce the correct binding energy, E cxp (^He) = —3.12 ± 0.02 MeV, with an 
adjustable parameter / = 0.8923. This overbinding problem is mainly attributed to 
the Brueckner rearrangement effect of the a-cluster, originating from the starting 
energy dependence of the bare two-nucleon interaction due to the addition of an 
extra A particle. 41 ) The odd-state AN interaction is assumed to be zero. All partial 
waves up to AMax = ^lMax = 6 for the aa and Aa pairs are included. The direct and 
exchange Coulomb kernel between the two a-clusters is introduced at the nucleon 
level with the cut-off radius, Rc = 14 fm. Table VIII shows the ground-state (0 + ) 
and the 2 + excitation energies of ^Be, predicted by the SB and the other various 
AN potentials used by Hiyama et a/. 42 ) In the present calculations using only the 
central force, the SB potential with the pure Serber character can reproduce the the 
ground-state and the excitation energies within 100 - 200 keV accuracy. 

4.4. The AAa system for A 6 A He 

We can use the Aa T-matrix, used in the aaA Faddeev calculation, to calculate 
the ground-state energy of XiEfe- The full coupled-channel T-matrices of fss2 and 
FSS with the strangeness S = — 2 and the isospin 1 = are employed for the AA 
RGM T-matrix. Table IX shows the ABaa values in MeV, predicted by various 
combinations of the AN and AA interactions. The results of a simple three-range 
Gaussian potential, Vaa (Hiyama), used in Ref. 42) are also shown. We find that 
this AA potential and our Faddeev calculation using the FSS T-matrix yield very 
similar results with the large ABaa values about 3.6 MeV, since the AA phases 
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Table IX. Comparison of ABaa values in MeV, predicted by various AA interactions and Van po- 
tentials. The AA potential Vyiyi(Hiyama) is the three-range Gaussian potential used in Ref. 42). 
FSS and fss2 use the AA RGM T-matrix in the free space. The heading eaa stands for the AA 
expectation value determined self-consistently, and Vaa{SB) the two-range Gaussian potential 
given in the text. In Vaa(SB) only the 5-wave is used, while in the others converged results 
with enough partial waves are given. 





Vaa (Hiyama) 


FSS 




fss2 




Vaa(SB) 


Van 


present 


Ref. 42) 


ABaa 


EAA 


ABaa 


eaa 


(S-wave) 


SB 


3.618 




3.650 


5.131 


1.411 


5.942 


1.902 


NS 


3.548 


3.59 


3.623 


5.159 


1.364 


5.951 


1.917 


ND 


3.181 


3.10 


3.231 


4.482 


1.288 


5.231 


1.636 


NF 


3.208 


3.22 


3.301 


4.626 


1.270 


5.410 


1.719 


JA 


3.370 


3.44 


3.468 


4.908 


1.306 


5.705 


1.833 


JB 


3.486 


3.56 


3.592 


5.150 


1.325 


5.956 


1.924 



shifts predicted by these interactions increase up to about 40°. The improved quark 
model fss2 yields ABaa = 1-41 MeV. (If we use the AA single-channel T-matrix, this 
number is reduced to 1.14 MeV.) In Table IX, results are also shown for Vm(SB), 
which is a two-range Gaussian potential generated from the fss2 1 S'o AA phase shift 
(see Fig. 11) in the full-channel calculation, using the supersymmetric inversion 
method. 39 ^ This potential is explicitly given by 

V^(SB) = -103.9 exp(-1.176 r 2 ) + 658.2 exp(-5.936 r 2 ) , 

where Vati(SB) in MeV and r in fm is the relative distance between two A's. We 
think that the 0.5 MeV difference between our fss2 result and the Vm(SB) result is 
probably because we neglected the full coupled-channel effects of the AAa channel 
to the SNa and EE a channels. We should also keep in mind that in all of these 
three-cluster calculations the Brueckner rearrangement effect of the a-cluster with 
the magnitude of about —1 MeV (repulsive) is very important. 27 ) It is also reported 
in Ref. 44) that the quark Pauli effect between the a cluster and the A hyperon 
yields a non-negligible repulsive contribution of 0.1 - 0.2 MeV for the A separation 
energy of ^He, even when a rather compact (3q) size of b = 0.6 fm is assumed as 
in our quark-model interactions. Taking all of these effects into consideration, we 
can conclude that the present results by fss2 are in good agreement with the recent 
experimental value, AB^ = 1.01 ± 0.20 MeV, deduced from the Nagara event. 26 ) 

§5. Summary 

From the advent of Yukawa's meson theory, a huge amount of efforts have 
been devoted to understand the fundamental nucleon-nucleon (NN) interaction 
and related hadronic interactions. The present view of these interactions are non- 
perturbative realization of inter-cluster interactions, governed by the fundamental 
theory of the strong interaction, quantum chromodynamics (QCD), in which the 
gluons are the field quanta exchanged between quarks. On this basis, meson "the- 
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ory" can be understood as an effective description of quark-gluon dynamics in the 
low-energy regime. The short-range part of the NN and hyperon-nucleon (YN) 
interactions are still veiled with unaccessible mechanism of quark confinement and 
multi-gluon effects. 

Here we have applied a naive constituent quark model to the study of baryon- 
baryon interactions, in which some of the essential features of QCD characteristics 
are explicitly incorporated in the non-relativistic framework. For example, the color 
degree of freedom of quarks is explicitly incorporated, and the full antisymmetriza- 
tion of quarks is carried out in the resonating-group (RGM) formalism. The gluon 
exchange effect is represented in the form of the quark-quark interaction, for which 
a color analogue of the Fermi-Breit (FB) interaction is used with an adjustable pa- 
rameter of the quark-gluon coupling constant as- The confinement potential is a 
phcnomenological r 2 -type potential, which does not contribute to the baryon-baryon 
interactions in the present framework. Since the meson-exchange effects are the non- 
perturbative aspect of QCD, these are described by the effective meson exchange 
potentials (EMEP) acting between quarks. After several improvements, the most 
recent model called fss2 has achieved very accurate description of the NN and YN 
interactions. x )~ 4 ) 

We have extended this (3q)-(3q) RGM study of the the AW and YN interactions 
to the strangeness S = —2, — 3 and —4 sectors without introducing any extra pa- 
rameters, and have clarified some characteristic features of the B$B$ inter actions. 5 ) 
The results seem to be reasonable, if we consider, i) the spin-flavor SUq symmetry, 
ii) the weak pion effect in the strangeness sector, and iii) the effect of the flavor 
symmetry breaking. These BgBg interactions are now used for the detailed study of 
the few-body systems, as well as baryonic matter problems, in various ways. Here we 
discussed applications of the NN, YN and YY interactions to the Faddeev calcula- 
tions of the three-nucleon bound state, the ANN-ENN system for the hypertriton, 
the aaA system for 9 Be, and the AAet system for J^He. We find that fss2 gives the 
triton binding energy large enough to compare with the experiment, without reduc- 
ing the deuteron D-state probability. 7 )' 29 ) The charge root-mean-square radii of 3 H 
and 3 He are also correctly reproduced. The application to the hypertriton calcula- 
tion shows that fss2 gives a reasonable result similar to the Nijmegen soft-core model 
NSC89, 35 ) except for an appreciable contributions of higher partial waves. 31 ) In the 
application to the aaA system, the aa RGM kernel with the three-range Minnesota 
force 38 ) and some appropriate AN force generated from the low-energy phase-shift 
behavior of fss2 yield the the ground-state and excitation energies of ^Be within 100 
- 200 keV accuracy. 37 ) The weak AA force AB%% = 1.01 ± 0.20 MeV deduced from 
the Nagara event 26 ) for ^He is reasonably reproduced by the Faddeev calculation 
of the AAa system, using the present Aa T-matrix and the full coupled-channel 
AA-EN-EE f -matrix of fss2. 

It is important to note that the newly developed three-cluster Faddeev formalism 
using two-cluster RGM kernels opens a way to solve few-baryon systems interact- 
ing by the quark-model baryon-baryon interactions without spoiling the essential 
features of the RGM kernels; i.e., the non- locality, the energy dependence and the 
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existence of the pairwise Pauli-forbidden state. It can also be used for the three- 
cluster systems involving a-clusters, like the ^Be and ^He systems. A nice point 
of this formalism is that the underlying NN and YN interactions are more directly 
related to the structure of the hypernuclei than the models assuming simple two- 
cluster potentials. In particular, we have found that the model fss2 yields a realistic 
description of many three-body systems including the three-nucleon bound state, the 
hypertriton, ^Be and ^He. 
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